Ab initio shallow acceptor levels in gallium nitride 
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Impurity levels and formation energies of acceptors in wurtzite GaN are pre- 
dicted ab initio. Beoa is found to be the shallow (thermal ionization energy ^ 
0.06 eV); Mgca and Znoa are mid-deep acceptors (0.23 eV and 0.33 eV respec- 
tively); Caca and Cdca are deep acceptors (~0.65 eV); SIn is a midgap trap with 
high formation energy; finally, contrary to recent claims, Cn is a deep acceptor 
(0.65 eV). Interstitials and heteroantisites are energetically not competitive with 
substitutional incorporation. 

1. Introduction 

Gallium nitride has recently emerged as a serious candidate for the fabrication 
of green to violet light emitting devices. Blue to yellow-green LED's0'i and, most 
recently, a GaN-based blue laser have been demonstrated! A basic problem that 
needs to be solved is that of reproducible and efficient p-tvpe doping of GaN. Mg was 
succesfully employed as p-dopant for device application£l only after elaborate post- 
growth treatments, and despite its large (~0.25 eV) ionization energy. The relevant 
basic mechanisms are still under scrutiny.cl 

There is thus an obvious interest in shallower acceptors for GaN, and first-principle 
theoretical investigations can reliably indicate suitable candidates to this end. In 
this paper we present a preliminary account of ab initio studies on impurities from 
Groups IIA, IIB, and IV in wurtzite GaN. Our main results are: Beca? with a thermal 
level at 0.06 eV above the valence band, is the shallowest acceptor reported so far 
in GaN; Mgca and Znca are found to be mid-deep acceptors (0.23 eV and 0.33 eV 
respectively) in good agreement with experiments; Cn, Caca, and Cdca are deep 
acceptors (around 0.65 eV); SIn is a midgap trap (1.97 eV) with high formation 
energy. Importantly, preliminary results on interstitials and heteroantisites indicate 
that these are energetically not competitive with substitutional incorporation. 

2. Method 

The calculation of acceptor levels and formation energies is non-trivial, especially 
in a computationally demanding system such as GaN. We calculate local-density- 
functional total energies and forces with a conjugate- gradient minimization scheme,i 
plane-waves, and ultrasoft pseudopotentials.i By the latter, we describe accurately 
all "difficult" electronic states (such as Ga and Zn 3c?, N 2p, Be 2s, Cd Ad, Mg and 
Ca semicore s and p, etc.) with a cutoff of 25 Ryd. The defects are simulated by 
single impurities in 32-atom wurtzite GaN supercells. Full geometry optimization is 
performed in all cases, using a newly developed method (details given elsewhere). 
The impurity formation energy in the charge state Q {Q is the number of electrons 



transferred from the defect to a reservoir of chemical potential fie) is given by 

EfiQ) = E*°*(Q) - n^G'^ - n^/i^ - + Qif^e + E,), (1) 



where E^°^{Q) is the defected supercell total energy in the specific charge state, /ig the 
electron chemical potential, Ey the top valence band energy, n*^^, n^, and are the 
number of Ga, N, and X atoms, and /i^*^, /i^, and //"^ their chemical potentials. These 
potentials must satisfy the equilibrium conditions with GaN and with appropriate 
solubility-hmiting X-N or X-Ga compounds (whose properties are also calculated ab 
initio). To study the case of maximum impurity solubility, we fix fix at the maximum 
value compatible with the solubility limit; in this case a single free potential remains 
(e.g. /zn), which is determined by the imposed growth conditions. 

The thermal ionization level e(0/— ) of a single acceptor is by definition the forma- 
tion energy difference of the two charge states Q=0 and Q=-l at /ie=0, or otherwise 
stated the electron chemical potential value at which the negatively charged state 
becomes favored over the neutral one. This effectively corresponds to the capture 
of an electron in the center with concurrent geometric relaxation. The optical level, 
accessible in luminescence experiments, is obtained summing the calculated Franck- 
Condon shift to the thermal energy. To refer the energies of different charge states to 
the same zero, we align the supercell local potentials in appropriate bulk-like regions. 
By bulk-like region we mean one where the macroscopic average^ of the local potential 
is flat and unaffected by the impurity. Such regions are chosen on the wurtzite c-axis, 
along which the impurity state decays faster (see below), and the inter-acceptor dis- 
tance is maximal (~ 20 bohr). Fig. 1 exemplifies the procedure for Zn: a plateau 
appears halfway between the impurities (placed at and 1), and it is there that we 
pick the potential values for the alignment. 

Ideally, one should reach bulk-like re- 
gions with the same potential value no mat- 
ter in which direction one moves away from 
the impurity; an upper bound of the align- 
ment uncertainty is thus the difference of 
the constant values along e.g. the a and 
c axis. This directional deviation is essen- 
tially all of our "formal" error bar in the 
calculated levels, as our calculations are 
otherwise well converged in terms of tech- 
nical ingredients. It amounts to ~ 0.05 eV 
for Zn (Fig. 1), and is generally of order 0.1 
to 0.2 eV at worst. 
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Fig. 1: Local potential macroscopic aver- 
ages on c-axis for bulk GaN (solid) and Zn° 
(dashed), and on a-axis for Zn" (dotted). 



3. Results and Discussion 

In Table 1, we summarize the results for the acceptors of groups II-A, II-B, and 
IV. A new candidate p-dopant emerges: Beca, with a thermal level at E^ + 0.06 



eV. Mg and Zn are progressively deeper (0.23 eV and 0.33 eV), while all others are 
definitely deep. In n-type conditions (/le ~ -Egap — 3.4 eV) the formation energies of 
all acceptors in their negative charge state (except SIn) are low or negative, which 
makes them efficient donor compensators. While in p-type and Ga-rich conditions 
compensation by nitrogen vacancy donors is known to be a problem, from calculations 

Table 1. Acceptors in wurtzite GaN. Formation energies E/ of the neutral charge state in N-rich 
conditions, thermal impurity levels, relaxation contribution E^, change in distance to neighbors in 
the a-plane (Ada) and along the c-axis (Adc) in percentage of the ideal bond length, solubihty hmits 
and their formation enthalpy. All energies are in eV. 
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Beca 


2.29 


0.06 


0.40 


—5.5 


-8.4 


BesNa 


6.70 


MgGa 


1.40 


0.23 


0.21 


+3.2 


+3.2 


MggNa 


4.96 
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0.33 
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+1.1 
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Zn3N2 


1.28 
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2.15 
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1.87 
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CasNa 
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Cdca 
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0.65 


1.40 


+9.2 


+8.8 


h-Cd 




Cn 


4.24 


0.65 


0.04 


+0.5 


-0.6 


d-C 




SIn 


5.49 


1.97 


3.36 


+12.0 


+11.8 


d-Si 






2.32 




0.07 


0.0 


+1.2 


w-GaN 


1.65 



of impurity concentrations accounting for the formation of singly-positive nitrogen 
vacancies (formation energies: 2.32 eV (N-rich) and 0.67 eV (Ga-rich) at fie = 0) 
and for overall neutrality,^ we find that acceptor-vacancy compensation in p-type 
conditions is essentially suppressed for all species in N-rich conditions. In p-type 
conditions, formation energies are low for Zn and Mg, somewhat higher for Be, high 
for Cn, very high for SIn. The carrier concentration at room temperature for the 
shallowest acceptor (Be) in N-rich conditions and at MOCVD growth temperatures, 
ranges from ~3xl0^^ cm~^ to ~4xl0^^ cm~'^ for a formation entropy of to 10 k^. 
The incorporation mechanism suggested^ for Mg may be helpful also in the case of 
Be: the calculated formation energy of the Be-H complex (limits: H2 and BesNa) 
is as low as 0.32 (0.87) eV/pair in N-rich (Ga-rich) conditions. The concentration 
of Be thus incorporated would grow by 3 orders of magnitude. Of course H is also 
incorporated, and must be detached from the complex to activate the acceptor: this 
requires about 1.5 eV for the Be-H complex, so that LEEBI or thermal annealing 
may be needed. 

Contrary to recent claims! on the Cn acceptor being shallow (~ 0.2 eV), we find 
it to be deep, with a thermal ionization energy of 0.65 eV. Due to its high formation 
energy in p-type conditions, any appreciable role of Cn as a shallow dopant, and in 
p-type conditions, is to be ruled out. On the other hand its formation energy drops to 
near zero in n-type conditions, so it may be present in high concentrations in n-GaN. 
The optically active level is 0.8 eV above the valence band. 

The electronic structure of the acceptors is unusual in several respects. First, the 



impurity states of the Ga-substituting species (including the mid-deep and shallow 
ones) are largely localized on the first-neighbor N shell. Second, the impurity states 
are of N p^^y-like character. This is a reminder of the impurity state's origin from 
the Fg top valence band doublet; as such, it is of course peculiar of wurtzite GaN. 
The state is generally localized in the a-plane, and decays rapidly along the c-axis. 
We display it for Mg in Fig. 2. The unusual shape of the impurity state affects the 
properties of H-acceptor complexes: e.g., H's preferential binding at the N-antibon- 
ding site in the a-plane rather than along the c-axis, correlates with the a-plane 
localization of the state accepting the extra electron. Third, in all cases, we find a 
sizable polarization of the ?>d states of second-neighbor Ga's. For Zn, the So? states 
are even shallower, and they mix into the impurity state. This may be the cause of 
the larger Zn binding energy in GaN as compared to e.g. GaAs. In the latter, the 
delocalized carrier bound to the shallow center sees the Ga and Zn cores as effectively 
identical: this may not hold in GaN, where the 3(i's hybridize with sp valence states. 

Finally, whether or not acceptors in GaN 
behave effective-mass like is an open issue. 
It may be speculated that the answer is in 
the affirmative: using an average hole mass 
m/i ~ 0.8 me the hydrogenic binding en- 
ergy is 0.4 eV (not far from the "isocoric" 
acceptor Zn) when the high-frequency di- 
electric constant £ = £oo = 5.4 is used. This 
is suggested by the fact that most accep- 
tor binding energies are much larger than 
typical GaN phonon energies.0 As for the 
shallower Be center. Be is expected to have 
repulsive central cell corrections, and its 
smaller binding energy may imply using a 
low-frequency dielectric constant £ = £:o~10 instead: in this case the binding energy 
is in fact ~ 0.1 eV. Of course the valence band structure and wave- vector-dependent 
screening must be described in detail before any conclusion can be reached. 
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Fig. 2: Number density of the impurity 
state for Mg. Note the Pxy-^i^e character 
and the locahzation on the first-neighbor 
shell. 
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